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ABSTRACT 

We report new observations of Zn II and Cr II absorption lines in 10 damped 
Lyman a systems (DLAs), mostly at redshift Zabs ~ 2.5. By combining these results 
with those from our earlier survey (Pettini et al. 1994) and other recent data, we 
construct a sample of 34 measurements (or upper limits) of the Zn abundance relative 
to hydrogen [Zn/H]; the sample includes more than one third of the total number of 
DLAs known. 

The plot of the abundance of Zn as a function of redshift reinforces the two main 
findings of our previous study. (1) Damped Lyman a systems are mostly metal-poor, 
at all redshifts sampled; the column density weighted mean for the whole data set 
is [Zn/H] = —1.13 lb 0.38 (on a logarithmic scale), or approximately 1/13 of solar. 
(2) There is a large spread, by up to two orders of magnitude, in the metallicities we 
measure at essentially the same redshifts. We propose that damped Lyman a systems 
are drawn from a varied population of galaxies of different morphological types and 
at different stages of chemical evolution, supporting the idea of a protracted epoch of 
galaxy formation. 

At redshifts z ^2 the typical metallicity of the damped Lyman a systems is in 
agreement with expectations based on the consumption of H I gas implied by the recent 
measurements of f^DLA by Storrie-Lombardi et al. (1996a), and with the metal ejection 
rates in the universe at these epochs deduced by Madau (1996) from the ultraviolet 
luminosities of high redshift galaxies revealed by deep imaging surveys. There are 
indications in our data for an increase in the mean metallicity of the damped Lyman a 
systems from z > 3 to «i 2, consistent with the rise in the comoving star formation 
rate indicated by the relative numbers of U and B drop-outs in the Hubble Deep Field. 
Although such comparisons are still tentative, it appears that these different avenues 
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for exploring the early evolution of galaxies give a broadly consistent picture. 

At redshifts z < 1.5 DLAs evidently do not exhibit the higher abundances expected 
from a simple closed-box model of global chemical evolution, although the number of 
measurements is still very small. We speculate that this may be due to an increasing 
contribution of low surface brightness galaxies to the cross-section for damped Lyman a 
absorption and to the increasing dust bias with decreasing redshift proposed by Fall 
and collaborators. However, more DLAs at intermediate redshifts need to be identified 
before the importance of these effects can be assessed quantitatively. 

The present sample is sufficiently large for a first attempt at constructing the 
metallicity distribution of damped Lyman a systems and comparing it with those 
of different stellar populations of the Milky Way. The DLA abundance histogram 
is both broader and peaks at lower metallicities that those of either thin or thick 
disk stars. At the time when our Galaxy's metal enrichment was at levels typical of 
DLAs, its kinematics were closer to those of the halo and bulge than a rotationally 
supported disk. This finding is at odds with the proposal that most DLAs are large 
disks with rotation velocities in excess of 200 km s~^, based on the asymmetric profiles 
of absorption lines recorded at high spectral resolution. Observations of the familiar 
optical emission lines from H II regions, which are within reach of near-infrared 
spectrographs on 8-10 m telescopes, may help resolve this discrepancy. 

Subject headings: cosmology: observations — galaxies: abundances — galaxies: 
evolution — quasars: absorption lines 
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1. INTRODUCTION 

In the last twelve months there has been a dramatic increase in our ability to identify normal 
galaxies at z ~ 3, study their stellar populations, and measure the rates of star formation and 
metal production in the universe over most of the Hubble time (Steidel et al. 1996; Madau et 
al. 1996). The most prominent features in the spectra of field galaxies at high-redshift (as is the 
case in the ultraviolet spectra of nearby star-forming galaxies) are strong interstellar lines which 
are similar, both qualitatively (in the range of ionization stages seen) and quantitatively (in the 
strengths of the absorption), to those in damped Lyman a systems; this similarity is consistent 
with the view that this class of QSO absorbers traces the material available for star formation at 
z ^ 2 (e.g. Wolfe 1995). The connection between normal galaxies and damped Lyman a systems 
(DLAs) is a particularly important one to make and clarifying several aspects of this connection 
remains a priority. The reason is simple: QSOs with known DLAs are typically more than 5 
magnitudes brighter than a L* galaxy at the same redshift. Consequently, we will inevitably 
continue to rely mostly on QSO absorption line spectroscopy for the study of physical conditions 
in the early stages of galaxy formation. 

Since 1990 (Pettini, Boksenberg, &; Hunstead 1990) we have been conducting a survey of 
metallicity and dust in DLAs taking advantage of the diagnostic value of weak transitions of Zn II 
and Cr II. As explained in that paper (see also the critical reappraisal of the technique in Pettini 
et al. 1997), [Zn/HjQ is a straightforward measure of the degree of metal enrichment analogous to 
the stellar [Fe/H], while [Cr/Zn] reflects the extent to which grain constituents are removed from 
the gas phase and thereby gives an indication of the dust-to-metals ratio. The major results of the 
survey were reported in Pettini et al. (1994). From the analysis of Zn and Cr abundances in 17 
DLAs, mostly at z ~ 2, we concluded that the typical metallicity of the universe at a look-back 
time of ~ 13 Gyr (Hq = 50 km s"^ Mpc"^; go = 0.01) was Zdla = I/IOZ©. We further found 

^We use the conventional notation where [X/Y] — log (X/Y)— log (X/Y)q 
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that there is a considerable range — by up to two orders of magnitude — in the degree of metal 
enrichment reached by different damped Lyman a galaxies at essentially the same epoch, and that 
even at these early stages of galaxy formation dust appears to be an important component of the 
interstellar medium, leading to the selective depletion of refractory elements from the gas. 

A natural next step is to extend the Zn and Cr abundance measurements over a wider range 
of redshifts than that considered by Pettini et al. (1994) with the ultimate aim of identifying 
the emergence of heavy elements and dust in galaxies and following their build-up with time. To 
this end we have continued our survey since 1994; the full sample now consists of 34 DLAs, more 
than one third of the total number known (Wolfe et al. 1995). In this paper we present the new 
data and consider the conclusions that can drawn from the whole set of measurements of [Zn/H]; 
preliminary reports have appeared in conference proceedings (e.g. Pettini et al. 1995a; Smith et 
al. 1996). Our findings on the abundance of dust from consideration of the [Cr/Zn] ratio in the 
same sample have been reported separately (Pettini et al. 1997). Recently, Lu et al. (1996) have 
addressed similar questions from measurements of [Fe/H] in 20 DLAs using high-resolution echelle 
spectra acquired with the Keck telescope. These authors reach conclusions which are in agreement 
with those presented here regarding the emergence of heavy elements at high redshifts, although 
the analysis of [Fe/H] is complicated by the fact that this ratio, unlike [Zn/H], depends on hoth 
the metallicity and dust content of the interstellar medium. 

Before proceeding it is useful to point out that in cases where DLAs from the present sample 
have been reobserved with HIRES on Keck (Wolfe et al. 1994; Wolfe 1995; Prochaska h Wolfe 
1997a), [Zn/H] has been found to be in good agreement with the values measured in our survey, 
which is based on 4-m telescope data (see §2 below). While the exceptional quality of the Keck 
observations has made possible several new aspects of this work, including the study of the relative 
abundances of a wide range of elements and the analysis of the kinematics of the absorbing gas, 
the basic survey of metallicity in DLAs can be carried out satisfactorily with 4-m class telescopes. 
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The main reason for this is the optically thin nature of the Zn II and Cr II lines in most DLAs 
proposed by Pettini et al. (1990) and confirmed by subsequent Keck spectra. 

2. OBSERVATIONS AND DATA REDUCTION 

The new data reported in this paper consist of observations of 10 DLAs in 9 QSOs 
obtained between March 1994 and February 1996 (an additional candidate DLA from the low 
dispersion survey by Storrie-Lombardi et al. (1996b) — at 2:abs=3.259 in the Zem=4.147 BAL QSO 
1144—073 — was shown not to be a damped system by our higher resolution observations of the 
Lyman a absorption line). In Table 1 we have collected relevant information for the 10 DLAs; 
the references listed in column (4) are the papers where the damped nature of the absorber was 
first identified. The absorption redshifts measured from associated metal lines in our blue and red 
spectra are listed in column (5); with 6 new DLAs at ^abs > 2.5, we have tripled the number of 
absorbers in this redshift regime compared with our earlier sample. 

The observations, reduction of the spectra and derivation of Zn and Cr abundances followed 
the procedures described in Pettini et al. (1994) and the interested reader is referred to that paper 
for a detailed treatment. Briefly, the observations were carried out mostly with the double-beam 
cassegrain spectrograph of the William Herschel telescope on La Palma, Canary Islands; additional 
red spectra were secured with the cassegrain spectrograph of the Anglo-Australian telescope 
at Siding Spring Observatory, Austraha. At Zabs > 2.5 the Zn II AA2025.483, 2062.005 and 
Cr II AA2055.596, 2061.575, 2065.501 multiplets are redshifted longwards of 7175 A, where the 
quantum efficiency of CCDs falls with increasing wavelength. Using EEV and Tektronix CCDs 
we generally found it necessary to integrate for longer than ^ 20 000 s (column 8 of Table 1) in 
order to achieve S/N between 9 and 46 (column 9). With a spectral resolution of 0.75-1.1 A 
FWHM (column 7), the corresponding Scr detection limits for the rest frame equivalent widths 



- 7- 



of unresolved Zn II and Cr II absorption lines range from Wo{3a) = 66 to 14 mA (column 10). 
The final "depth" of the survey — that is the lowest metallicity that can be measured — depends 
on the combination of Wo{3a) and the neutral hydrogen column density A''(H''). Since the 
values of Ar(H°) in the new DLAs observed span one order of magnitude (see §3 below), it is the 
sight-lines with the largest column densities of gas which provide the most stringent limits on 
metal abundances. Accordingly, we have tended to select DLAs for the present survey primarily 
on the basis of the value of Ar(H°). 

In Figure 1 we have reproduced portions of the QSO spectra encompassing the regions where 
the Zn II and Cr II lines are expected in the 10 DLAs in Table 1. As can be seen from the figure, 
the absorption lines sought are detected in approximately half of the cases. Table 2 lists redshifts 
and rest-frame equivalent widths for the detections; in the other cases the 3a" limits given in 
column (10) of Table 1 apply. 

With the double-beam spectrograph on the WHT we were able to record portions of the blue 
spectrum of each QSO, centred on the damped Lyman a line, simultaneously with the red arm 
observations aimed at the Zn II and Cr II lines. The blue detector was either the Image Photon 
Counting System or a thinned Tektronix CCD; exposure times were the same as those given in 
column (8) of Table 1. A 600 grooves/mm grating was used to record a 800 A wide interval of the 
spectrum with a resolution of ~ 1.5 A FWHM. This configuration was chosen in preference to the 
higher resolving power achievable with a 1200 grooves/mm grating because a good definition of 
the QSO continuum is a key factor in determining the accuracy with which N(il^) can be deduced 
from the profile of the damping wings of the Lyman a absorption line. Normalised portions of 
the blue spectra are shown in Figure 2 together with our fits to the damped Lyman a lines. The 
theoretical damping profiles are centred at the redshifts of the O I A1302.1685 lines which are 
encompassed by our blue data. 
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3. ZINC AND CHROMIUM ABUNDANCES 

The main results of our survey are collected in Table 3 which includes the 10 DLAs in Table 
1 and 7 additional systems for which data have been published since our earlier study (Pettini et 
al. 1994). Values of the neutral hydrogen column density A'^(H'') are listed in column (3) of Table 
4; the typical accuracy of these measurements, including the uncertainty in the placement of the 
continuum, is ±20%. A/"(H°) is likely to account for most of the neutral gas in each DLA given the 
low molecular fractions which apply to these absorbers at high redshifts (Levshakov et al. 1992; 
Ge & Bechtold et al. 1997; Cirkovic et al. 1997). 

Columns (3) and (6) of Table 3 give the column densities of Zn"*" and Cr"*" respectively, 
deduced from the measured equivalent widths (or upper limits) assuming no line saturation. 
That this is generally the case is indicated by: (1) the weakness of the absorption lines; (2) the 
equivalent width ratios of lines within each multiplct which, when measurable, are usually close to 
the ratios of the corresponding /-values (Bergeson & Lawler 1993); and (3) the resolved absorption 
profiles recorded with HIRES on Keck for many DLA systems, including some in common with the 
present survey (Lu et al. 1996; Prochaska & Wolfe 1997a). There are of course exceptions, such as 
the 

-^abs — 2.5842 system in Q1209+093 — see the discussion at §3.12 below. The important point, 
however, is that it is usually possible with the signal-to-noise ratio and resolution of our data to 
assess the degree of saturation of the Zn II and Cr II lines. 

Column (4) lists the ratios A?^(Zn+)/A/"(H°) derived by dividing the entries in column (3) by 
those in column (3) of Table 4; comparison with the solar abundance of Zn, log (Zn/H)0 = —7.35 
(Anders & Grevesse 1989), then leads to underabundances of Zn by the factors given in column 
(5). The corresponding values for Cr (log (Cr/H)Q = —6.32) are given in column (8) and column 
(9) lists the ratio A?^(Cr+)/7V(Zn+) in cases where it could be determined. 

In taking the ratios N{Zn+)/N{li^) and N{Cv+)/N{'H.^) as measures of (Zn/H) and (Cr/H), 
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we implicitly assume that there is little contribution to the observed Zn II and Cr II absorption 
from ionised gas (which would not produce Lyman a absorption). This is likely to be the case 
given the large column densities of H I, and indeed there are no indications to the contrary in our 
data. In particular, we found no significant differences in redshift between the Zn II and Cr II 
lines, when detected, and O I A1302.1685 which arises only in H I regions. Should this assumption 
be shown to be incorrect, however, the values of [Zn/H] and [Cr/H] deduced here would then be 
upper limits to the true abundances. 

We now comment briefly on each DLA in Table 3. 

3.1. QOOOO-263; Zabs = 3.3901 

Our observations of this DLA, the highest redshift absorber in the survey, have been described 
in Pettini et al. (1995a). While Zn II A2025.483 remains undetected, despite the sensitive limit 
reached in a total exposure time of 58 200 s, we do record weak Cr II absorptions at the 4a 
(A2055.596) and 3a (A2061.575) significance levels. Cr II A2055.596 is expected to be stronger 
than Zn II A2025.483 if the fraction of Cr locked up in dust grains is less than about 50%. With 
AT (HO) = (2.5 ± 0.5) X 10^1 cm'^ (Savaglio, D'Odorico, & Moller 1994), this is one of the highest 
column density systems in our sample. We conclude that the abundance of Zn is less than 1/80 of 
the solar value; this estimate is ~ 5 times more sensitive than the previous limit (Savaglio et al. 
1994). The abundance of Cr, [Cr/H] ~ —2.2 it 0.1, is similar to those of other elements measured 
by Molaro et al. (1996) and Lu et al. (1996), making this DLA one of the most metal-poor in our 
sample. 
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3.2. Q0056+014; Zabs = 2.7771 

This QSO is from the Large Bright Quasar Survey by Chaffee et al. (1991). We deduce 
A^(H*') = (1.3 lb 0.2) X 10^^ cin~^ from fitting the core of the damped Lyman a Hne, in reasonable 
agreement with the value log A'^(H I) = 21.0 reported by Wolfe et al. (1995). 

As can be seen from Figure 1, the Zn II and Cr II absorption lines are broad and shallow 
in this DLA, spanning ~ 200 km s~^. The stronger member of the Zn II doublet, A2025.483, 
falls within the atmospheric A band. Plotting the four absorption lines labelled in Figure 1 on 
the same velocity scale suggests that most of feature "1" is not due to Zn II A2025.483, but 
rather to poorly corrected telluric absorption. From the equivalent widths of Cr II A2055.596 
and A2065.501 (features 2 and 4 in Figure 1), which are consistent with the optically thin ratio 
of 2:1, we deduce a weighted mean iV(Cr+) = (2.8 ± 0.4) X 10^^ cm^^ This column density of 
Cr+ produces an equivalent width Wq =(82 ib 12) mA for Cr II A2061.575; since we measure 
Wo =(117 ± 16) mA for feature 3, which is a blend of Cr IIA2061.575 and Zn IIA2062.005, 
we conclude that Wq =(35 ± 20) mA for the latter. This in turn corresponds to Ar(Zn+) = 
(3.5 lb 2)xl0^^ cm~^. Thus both Zn and Cr appear to be ^ 20 times less abundant than in the 
Sun. 

Our red spectrum also shows several Fe II lines from an absorption system at Zabs = 2.3044, 
including: Fe II A2344.214 (visible in Figure 1 at Aobs = 7748.46A) with Wq =(470 ± 12) mA; 
Fe II A2367.5905 with Wq =(64 ± 6) mA; Fe II A2374.4612 with Wq =(220 ± 14) mA; and 
Fe II A2382.765 with Wq =(640 ± 12) mA. 
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3.3. Q0201+365; Zabs = 2.462 

Keck observations of this DLA have been pubhshed recently by Prochaska & Wolfe (1996) 
who deduced relatively high abundances of Zn and Cr, respectively ^1/2 and ~ 1/8 of solar. 
Evidently, even at redshifts as high as 2.5 some galaxies had already undergone significant chemical 
evolution and enriched their interstellar media in heavy elements to levels comparable with that 
of the Milky Way today. 

3.4. Q0302-223; Zabs = 1 0093 

Lanzetta, Wolfe, & Turnshek (1995) proposed this as a candidate DLA system on the basis of 
low-resolution lUE data; a subsequent UV spectrum secured with the Faint Object Spectrograph 
on the Hubble Space Telescope confirmed that N(H.^) = (2.15 ± 0.35) x 10^° cm~^ (Pettini & Bowen 
1997). Recent WHT observations of Zn II and Cr II lines by Pettini &; Bowen (1997) have shown 
the abundances to be 1/3 and 1/8 of solar respectively. After subtraction of the QSO radial profile 
from HST WFFC2 images of the field, Le Brun et al. (1997) identified two galaxies which may be 
producing the absorption; at 2; = 1.009 they would have luminosities L « 0.2L* and ^ L* and 
distances of 12 and 27 h^Q kpc respectively from the QSO sight-line. 

3.5. Q0454-F039; Zabs = 0.8596 

The abundances of Zn and Cr reported by Steidel et al. (1995a) correspond to [Zn/H] = 
—0.83 lb 0.08 and [Cr/H] = —1.01 it 0.05 if the experimentally measured /-values of the Zn II and 

Cr II multiplets (Bergcson & Lawlcr 1993) are adopted for consistency with the rest of the present 
study. Deep images of the QSO field both from the ground (Steidel et al. 1995a) and with HST 
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(Le Brun et al. 1997) suggest that the absorber is a compact galaxy with L 0.25L* (go = 0.05) 
at a projected distance of 8 h'^Q kpc from the QSO. 

3.6. Q0836+113; Zabs = 2.4651 

This is the faintest QSO in our survey (Hunstead, Pettini, & Fletcher 1990) and the S/N of 
the red spectrum remains modest despite the considerable investment in exposure time (Table 
1). Combined with the relatively low H I column density of (3.8 ± 0.4)xl0^'^ cm~^, the 3cr upper 
limits to the Zn II and Cr II lines place limits on the abundances of Zn and Cr which are less 
stringent than in most other DLAs considered: [Zn/H] < —0.8 and [Cr/H] < —1.2. 

The blue spectrum shown in Figure 2 was recorded with the IPCS on the WHT in 
March 1994. Note that, of all the damped Lyman a lines reproduced in Figure 2, this is 
the only instance where there appears to be weak emission in the core of the absorption 
line. The line flux, (2 it 0.7)xl0~^^ erg cm~^, agrees within the errors with the value of 
(2.9 lb 0.7)xl0~^^ erg cm~^ reported by Hunstead et al. (1990) from independent data 
obtained in April 1987 with a different IPCS detector on the AAT. The two sets of observations 
were obtained with the same slit width (1.2 arcsec) and at the same position angle on the sky (150 
degrees) . 

3.7. Q0841+129; Zabs = 2.3745, 2.4764 

The spectrum of this bright {V 2± 17), high redshift {z 2± 2.5, estimated from the onset of 
the Lyman a forest) BL Lac object discovered by C. Hazard (private communication) shows two 
DLAs (see Figure 2), making it a highly suitable target for follow-up high resolution observations. 
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As can be seen from Figure 1, in the lower redshift system we detect features 2 and 3; the 
strength of the latter indicates a significant contribution from Zn II A2062.005 to the blend. 
Following a procedure similar to that described for Q0056+014 at §3.1 above, we deduce A'^(Cr+) 
= (9.5 ± 2) X 10^2 cm-2 from the equivalent widths of Cr II A2055.596 and A2065.501. This in 
turn leads us to estimate that approximately half of the equivalent width of feature 3 is due to 
Zn II A2062.005 with Wq =(24 ± 9) mA. Together with the 3cr upper limit Wo{2025)< 26 ml for 
the stronger member of the doublet, this then implies A'^(Zn+) = (1.8 =b 0.5) x 10^^ cm~^. 

Thus we find that Zn and Cr at Zabs = 2.3745 are underabundant by factors of 23 and 45 
respectively, relative to solar values. Similar, or lower, abundances apply to the Zabs = 2.4764 
DLA, given the lack of detectable Zn II and Cr II lines (see Table 3). 

3.8. Q0913+072; Zabs = 2.6183 

The signal-to-noise ratios of our spectra of this bright QSO are among the highest in the 
survey — see Table 1 and Figures 1 and 2. The column density of neutral hydrogen is however 
comparatively low, A'^(H'') = (2.3 ± 0.4) xlO^'^ cm~^. The lack of Zn II and Cr II absorption even 
at S/N = 46 implies underabundances by factors of more than 14 and 32 respectively. 

3.9. Q0935+417; ^abs = 1-3726 

Lanzetta et al. (1995) estimated iV(H°)~ 2 x 10^° cm'^ for this candidate DLA 
from low resolution lUE data; a subsequent HST FOS spectrum confirmed that N{lP) = 
(2.5 ± 0.5) X 10^° cm~^ (Lanzetta Sz Meyer 1996, private communication). With this value of the 
hydrogen column density, the observations by Meyer, Lanzetta, & Wolfe (1995) imply [Zn/H] = 
-0.80 and [Cr/H] = -0.90 . 
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3.10. Q1104-180; Zabs = 1-6616 

Smettc ct al. (1995) identified this DLA in the spectrum of the brigliter (B = 16.7) component 
of this gravitationahy lensed QSO pair. Prom AAT observations obtained with an instrumental 
setup similar to that used in our survey, these authors estimated iV(H°)= 6 x 10^° cm~^. 
They also reported detections of Zn II and Cr II absorption lines with equivalent widths 
W^o (2025.483)= (75 ± 20) mA and W^o(2055.596)= (57 ± 20) mA respectively. If the lines are 
unsaturated [Zn/H] = -0.80 and [Cr/H] = -1.30. 

3.11. Q1151+068; ^abs = 1-7736 

Even though the damped Lyman a line falls in the crowded near-UV spectrum of this 
Zem = 2.762 QSO (see Figure 2), our estimate Ar(H°) = (2.0 ± 0.5)xl02i cm-^ is in very good 
agreement with log A^(H I)= 21.3 published by Turnshek et al. (1989). The ratios of equivalent 
widths within the Zn II and Cr II multiplets strongly suggest that the lines are optically thin; Zn 
and Cr are both underabundant by a factor 40. 

Our red spectrum, which covers the region 5500 - 5900 A, shows three C IV AA1548, 1550 
doublets at Zghs = 2.5629, 2.7069 and 2.7551 respectively. 

3.12. Q1209+093; Zabs = 2.5843 

This is another high column density DLA; we measure N(H.^) = (2.0 it 0.5) x 10^^ cm~^ 
which compares well with log N{}i I)= 21.4 reported by Lu et al. (1993). The Zn II and Cr II 
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lines are the strongest encountered in the entire survey of 34 DLAs (see Table 2). Fitting the 
absorption profiles requires 6 = 50 km (as usual, b = V2a where a is the one-dimensional 
velocity dispersion along the line of sight), indicating that several velocity components most 
likely contribute to the absorption. Similarly, Lu et al. found b = 122 km s~^ from fitting a 
single curve-of-growth to the strongest UV absorption lines. Some of the components may well be 
saturated; so we quote our best estimates of A''(Zn"'") and A''(Cr"'") as lower limits. We conclude 
that Zn is more abundant than 1/9 solar and Cr more abundant than 1/27 solar. Higher resolution 
observations are required to establish how close to these limits the true values are. 

3.13. Q1328+307; Zabs = 0.6922 

We have included here the measurements of Zn and Cr abundances in the spectrum of 3C 286 
reported by Meyer Sz York (1992), after appropriate rescaling with the /-values of Bergeson & 
Lawler (1993). Although the measurement was discussed in Pettini et al. (1994), this intermediate 
redshift DLA was not explicitly included in that sample which consisted exclusively of DLAs 
at Zabs > 1-7. CCD images of the QSO field obtained with ground-based telescopes (Steidel et 
al. 1994) and with HST (Le Brun et al. 1997) show a large (« 10 — 30 h^Q kpc), low surface 
brightness galaxy which has been proposed as the absorber. 

3.14. Q1337-M13; Zabs = 2.7957 

Our measured iV(H°) of (8 ± 2)xl0^° cm ^ agrees very well with log A^(H I)= 20.9 reported 
by Turnshek et al. (1989). When we observed this QSO, in March 1994, we found it to be 
considerably fainter than the magnitude V 2± 18.2 estimated by Hazard et al. (1986) from POSS 
plates. Although the S/N achieved is the lowest in the survey (see Table 1 and Figure 1), it is 



still sufficient to establish that the abundances of Zn and Cr are less than 1/10 and 1/23 of solar 
respectively. 



3.15. Q1946+769; z^),^ = 2.8443 

This Zem = 3.051 QSO, intrinsically one of the most luminous known, is sufficiently bright 
to have been studied extensively at echelle resolutions and high S/N with 4-m telescopes (Fan & 
Tytler 1994; Lu et al. 1995; Tripp, Lu, & Savage 1996). However, the hydrogen column density in 
the Zabs = 2.8443 DLA is relatively low, iV(H°) = (2± 0.5)xl020 cm-^ (Lu et al.). Consequently, 
the upper limits [Zn/H] < —0.82 and [Cr/H] < —1.00 deduced by these authors are rather 
uninformative given that the true metallicity is ~ 30 times lower ([Fe/H]= —2.44 ± 0.13). 



3.16. Q2239-386; Zabs = 3.2810 

This QSO is faint and the absorber is at high redshift; the combination of these two 
factors resulted in the longest integration time in the survey (see Table 1). Adopting 
N{'H.^)= 5.8 X lO^'^ cm~^ measured by Lu & Wolfe (1994), we deduce Zn and Cr underabundances 
by factors of more than 11 and 13 respectively. 

The Cr measurement is based on the weakest member of the triplet, Cr II A2065.501; 
A2061.575 is affected by a strong sky emission line and A2055.596, which at z^hs = 3.2810 is 
redshifted to Aobs = 8802.82 A, falls very close to Mn II A2606.462 at Zabs = 2.3777, the redshift of 
a second DLA along this line of sight (Lu &; Wolfe 1994). Based on the strengths of the other two 
members of the Mn II triplet, A2576.877 at Aobs = 8703.55 A and A2594.499 at Aobs = 8763.97A, 
the feature labelled 2 in the last panel of Figure 1 is mostly Mn II A2606.462 . The two strong 
absorption lines also visible in this figure are Fe II AA2586.6500, 2600.1729 at Zabs = 2.3777. 
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4. DISCUSSION 

Our total sample, which consists of measurements (or upper limits) of [Zn/H] in 34 DLAs 
over the redshift range Zabs = 0.6922 — 3.3901, is constructed by combining data for the 17 DLAs 
in Table 3 with those for the 15 DLAs in Table 3 of Pettini et al. (1994) and with the further 
addition of two DLAs in Q0528— 250 (Meyer et al. 1989) which were included in the sample 
considered by Pettini et al. (1994) but not listed in their Table 3. All the points in Figure 3 are 
based on the /-values of the Zn II doublet by Bergeson & Lawler (1993) and the meteoritic solar 
abundance of Zn from the compilation by Anders Sz Grevesse (1989)0. 

We now consider what implications can be drawn from this extensive survey on the chemical 
evolution of the neutral content of the universe and on the relationship of damped Lyman a 
systems to present-day spiral galaxies. 

4.1. Chemical Evolution of Damped Lyman a Systems 

Figure 3 shows the abundance of Zn as a function of redshift. The enlarged sample confirms 
the two main conclusions reached by Pettini et al. (1994): 

(1) Damped Lyman a systems, at all redshifts probed, are generally metal-poor and presumably 
arise in galaxies at early stages of chemical evolution. 

(2) There appears to be a large range in the values of metallicity reached by different galaxies 
at the same redshift, pointing to a protracted 'epoch of galaxy formation' and to the fact that 

^This set of atomic parameters and solar abundance introduces a correction of —0.148 to the values of [Zn/H] 
published in Pettini et al. (1994) 
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chemical enrichment probably proceeded at different rates in different DLA galaxies. 

Wfiile we find gas with near-solar metallicities at redshifts as high as 2; ~ 2.5, there are also 
examples of galaxies with abundances less than 1/10 solar at a time when the disk of the Milky 
Way differed little from its present-day composition. At redshifts z ~ 2 — 2.5 the full range of metal 
abundances spans about two orders of magnitude. Although for metallicities Zdla ^ 1/50Zq the 
Zn II lines become vanishingly small and only upper limits to the abundance of Zn can be deduced, 
we do know from echelle spectroscopy of more abundant astrophysical elements that values of 
Zdla ^ 1/100Zq are not uncommon at Zabs = 2 — 3 (see Figure 1 of Pettini et al. 1995a). 

These two results are considered quantitatively in Table 5 where in the last column we list, 
for various subsets of our sample, the column density-weighted mean abundance of Zn 

[(Zn/HoLA)] = log((Zn/H)DLA) " log (Zn/H)^, (1) 

where 

((Zn/H)i,LA) = — , (2) 

1=1 

and a', the standard deviation from the column density weighted mean, defined as 

n 

in - 1) a'' = i:([ZVH], - [(Zn/Hj^LA)])' (3) 

i=l 

The summations in equations (2) and (3) are over the n DLA systems considered in each 
subsample. 

Under the working assumption that DLAs account for most of the material available for 
star formation at high redshift, the quantity [(Zu/Hj^la)] is a measure of the degree of metal 
enrichment reached by the universe at a given epoch. This is a general statement which follows 
from the column density distribution of Lyman a systems (Lanzetta et al. 1995) and which holds 
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irrespectively of the precise nature of the damped absorbers (disks, spheroids, gas clouds yet to 
collapse into galaxies, etc.), provided that there are no significant biases in the samples of DLAs 
from which our observations are drawn (Fall 1996). 

The values of [(Zn/Hj^^A)] ™ Table 5 are strictly upper limits (with the exception of subsample 
number 1), since the averages include systems for which only upper limits to the abundance of 
Zn are available. However, we expect the corrections to be small because the systems where 
the Zn II doublet is below our detection limits are preferentially those with the lowest values of 

n 

hydrogen column density A'^(H''). Specifically, the fractions of N{ii )^ contributed by DLAs 

1=1 

with undetected Zn II lines are 28% for the full sample, and 16%, 16%, and 37% respectively for 
subsamples 2, 3 and 4. To show that including the upper limits as detections has only a modest 
effect on the mean values of metallicity deduced, we have recalculated [(Zn/Hj-)L^)] for the full 
sample twice, substituting 2(j and la limits respectively in place of the 3(7 limits used in Table 
5 (it could indeed be argued that 3cr limits for the entire ensemble on Zn non-detections is an 
overly conservative approach). In this case, [(Zn/Hj-jLA)] decreases from —1.13 it 0.38 (the value 
listed in Table 5) to —1.16 ± 0.40 and —1.20 ± 0.48 respectively. On the other hand, all three 
measurements in subsample 5 (^^abs = 3.0 — 3.5) are upper limits and accordingly we quote the 
value of [(Zn/Hj3L^)] in this redshift interval as an upper limit. 

For the full sample of 34 DLAs in the range Zabs = 0.7-3.4 we find [(Zn/H^LA)] = -1.13±0.38. 
This is the same value as obtained by Pettini et al. (1994) when account is taken of the different 
/-values and solar abundance scale used in our earlier study. For comparison, [Zn/H]gas = —0.19 
along unreddened sight-lines in the solar vicinity (Roth &; Blades 1995; Sembach et al. 1995 — both 
analyses used the same /-values and solar scale as here) . If the interstellar medium (gas-|-dust) 
near the Sun has the same composition as the Sun, this would imply that approximately 35% 
of Zn is in solid form. On the other hand, Pettini et al. (1997) found that for the present 
sample of DLAs the typical dust-to-metals ratio is approximately half that of the Galactic ISM. 
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If we assume, therefore, that on average 83% of Zn in DLAs is in the gas phase, we obtain 
[(Zn/Hj5LA)] = -1-13 ± 0.38 + log(l/0.83) = -1.004 ± 0.38, and conclude that the column density 
weighted abundance of Zn in DLAs is 1/11 of that of the Milky Way ISM today. 

One of the motivations of the present work was to determine the redshift evolution of the 
metallicity of DLAs and thereby trace the increase of heavy elements in the universe from the 
epoch of galaxy formation to the present time. Prom Figure 4, where our measures of [(Zn/H^LA)] 
from Table 5 are plotted versus redshift, it can be seen that any such evolution is only mild in the 
present sample. Between z = 3 and 1.5, to which 80% of the sample refers, there appears to be 
little change from the typical [(Zu/Hj-jla)] = —1.13. This is less surprising, however, when one 
considers that this redshift interval spans a period of only 3 Gyr from 14.3 to 11.4 Gyr ago 
(Ho = 50 km Mpc~^; qq = 0.01) and that at these epochs evidently there was a large spread 
in the chemical enrichment of different DLA galaxies. 

On the other hand, the upper limit [(Zn/Hj^L^)] < —1.39 at z > 3 is lower than the means 
in the other redshift bins, providing tentative evidence for a rapid build-up of elements with time 
at this epoch. This suggestion is strengthened by the data of Lu et al. (1996) who found that 
[Fe/H] ^ —2 in three additional DLAs at 2; > 3 . (The correction to [Fe/H] for the fraction of Fe 
in solid form is likely to be small — probably less than a factor of two — at such low metallicities; 
Pettini et al. 1997). The lowest metallicities measured in DLAs, Zdla — —2.5, are comparable 
to those thought to apply to the ionised intergalactic medium producing the Lyman a forest at 
redshifts z = 2 — 3.5 (Hellsten et al. 1997), although the large ionization corrections involved make 
estimates of Ziqm considerably more uncertain than ^dla- It is tempting, therefore, to interpret 
the rapid increase in metal abundances at 2; < 3 as an indication of the onset of star-formation 
in galaxies and to speculate that Z ~ —2.5 may be an approximate 'base' level of metallicity on 
which galactic chemical evolution subsequently builds. 

The recently realised ability to image high redshift galaxies directly in their ultraviolet stellar 
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continua has led to the first attempts to sketch the global history of star formation over ~ 80% 
of the age of the universe (Madau et al. 1996 and references therein). Determinations of the 
volume-averaged star formation rate (SFR) from the so-called B and U drop-outs (galaxies with 
the Lyman limit in the B and U bands respectively) in ground-based surveys (Steidel, Pettini, 
Sz Hamilton 1995b) and in the Hubble Deep Field (Madau 1996) do indeed suggest an increase 
in the SFR from z ~ 4 to z ~ 2.75 . As discussed by Madau et al., it is possible to convert the 
integrated UV luminosity density into a metal ejection rate p'z per comoving volume at redshift 
z. Since the massive stars which are the main contributors to the far-UV continuum are also the 
major producers of heavy elements (at least those released into the ISM by Type II supernovae), 
the conversion does not depend sensitively on the shape of the IMF in these primordial galaxies. 
Rather, the principal sources of uncertainty arise from the cosmology assumed and from the 
amount of dust extinction suffered by the UV continuum. 

Bearing in mind these uncertainties, it is of great interest to compare the values of 
Zdla deduced here with the metallicities which may be expected on the basis of Madau's 
metal ejection rate. Integrating p'z in Figure 3 of Madau (1996) from z = 5.5 to the present 
time yields a total density of metals pz{z = 0) 6.2 x 10^ Mq Mpc~^. This corresponds 
to an approximately solar metallicity if the present day density of baryons in galaxies is 
p^z = 0) fs (2.7 ± 0.4) X 10® M0 Mpc'^, or fi* « 4 x lO'^ {Ho = 50 km s"^ Mpc"^; Madau et 
al. 1996). In a closed box model, assuming that Cl^{z = 0) ^dla{z = 4) (Storrie-Lombardi, 
McMahon, & Irwin 1996a) we can take 



provided the gas consumption into stars from z = 5.5 to z is low and ^gas{z) S> 0.^{z). The 
redshift evolution of fioLA (Storrie-Lombardi et al. 1996a) suggests that this may well be the case 
up to z ~ 2 (as we proposed in Pettini et al. 1994). 




(4) 
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The broken line in Figure 4 shows the increase of Z{z)/Zq with decreasing redshift calculated 
from equation (4) and Madau's (1996) estimates of pz- Evidently, there is rough agreement 
between the predicted and observed values of ^dla- Given the current uncertainties, we consider 
it premature to read too much into this comparison. For example, Madau's p'z refers primarily to 
oxygen and the a-elements which presumably are more abundant than zinc and iron by a factor 
of 2 — 3 at these low metallicities (Edvardsson et al. 1993; Carney 1996). On the other hand, 
the broken line in Figure 4 may well underestimate the metal production rate by similar factors 
if star-forming galaxies at high redshift are reddened by small amounts of dust, corresponding to 
E(B — V) « 0.1, as suggested by the observed slopes of the UV continua (Steidel et al. 1996). 

Nevertheless, taken at face value. Figure 4 does seem to indicate that in the DLAs we see 
roughly the same level of metal enrichment as expected from direct observations of star-forming 
galaxies at these redshifts. More complex galactic chemical evolution models which use as a 
starting point the gas consumption indicated by the redshift evolution of JIdla (Pei & Fall 1995; 
Fall 1996) also reproduce the degree of metal enrichment of DLAs and the comoving rate of 
star formation at z ^2 . Thus it appears that, to a first approximation at least, these three 
independent avenues to exploring the epoch of galaxy formation — the consumption of neutral gas, 
the metal abundance of the absorbers, and the UV luminosity of high-redshift galaxies — give a 
broadly consistent picture of the early evolution of galaxies. 



4.1.1. Abundances at z < 1.5 

The situation is less clear at lower redshifts. Only four measurements make up subsample 
1 in Table 5, even though this bin spans a larger interval of time than all the other subsets put 
together — 5 Gyr from 11.4 to 6.3 Gyr ago (again for Hq = 50 km s~^ Mpc~^; qo = 0.01). 
Evidently [(Zu/Hj^la)] ~ —0.98 it 0.33 is below an extrapolation of Madau's curve in Figure 4. 
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However, it is difficult to assess liow firm tfiis conclusion is, given that 65% of E N{ir)i for 

i=l 

subsample 1 is due to the Zabs = 0.6922 absorber in 3C 286 which appears to be a large, low surface 
brightness galaxy (see §3.13 above). Possibly such galaxies, whose low metallicities at the present 
time are thought to be the result of low star formation efficiencies (McGaugh 1994; Padoan, 
Jimenez, &: Antonuccio-Dclogu 1997), come to dominate the cross-section for DLA absorption at 
z ^ 1, if by this epoch most high surface brightness galaxies have already processed a significant 
fraction of their gas into stars. Furthermore, the build-up of dust which goes hand-in-hand with 
the production of metals is likely to introduce an increasing bias (with decreasing redshift) against 
galaxies in advanced stages of chemical evolution, since existing samples of damped Lyman a 
systems are mostly drawn from magnitude limited optical QSO surveys (Fall &; Pei 1993; Pei &; 
Fall 1995). 

At Zabs < 1-5 imaging of DLA absorbers is within current observational capabilities. Although 
positive identifications based on spectroscopic redshifts have not yet been achieved, the candidates 
which have been proposed suggest a very diverse population of galaxies. While in some cases the 
absorbers could be on evolutionary paths similar to that of the Milky Way, the Zabs = 1-0093 DLA 
in Q0302— 223 being a good example (Pettini & Bowen 1997), there are also several instances 
where galaxies of low luminosity {Lb ^ O.IL*) or of low surface brightness are indicated (Steidel 
et al. 1994, 1995, 1997; Le Brun et al. 1997). 

Thus both effects considered above — a shift of the DLA population away from 'normal' L* 
galaxies and an increasing dust bias — may contribute to the finding that .^dla does not increase 
significantly at Zabs < 1-5 in Figures 3 and 4, contrary to simple expectations in a closed-box model 
of chemical evolution. However, it will not really be possible to proceed further, and quantify the 
relative importance of these two effects, without a larger sample of DLAs at intermediate redshifts. 
Identifying such a sample remains an urgent priority. 
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4.2. Comparison with Stellar Populations of the Milky Way 

Damped Lyman a systems are commonly thought of as the high redshift counterparts of 
present-day galactic disks, although we and others (Pettini et al. 1990; York 1988) have often 
made the point that high column densities of neutral gas are not the prerogative of disk galaxies 
alone. The sample of [Zn/H] measurements now available is sufficiently large to allow a comparison 
to be made of the distribution of metallicities in DLAs with those of different stellar populations 
in the Milky Way. In the solar cylinder, stars in the halo, thick disk, and thin disk have distinct 
dynamical and abundance properties, although the distributions overlap in either parameter 
taken separately. It is the combination of chemical abundance and kinematic data that studies of 
Galactic evolution have focussed on; here we attempt to use this information to throw light on the 
nature of DLA galaxies. 

Our measurements of [Zn/H] from Figure 3 have been plotted in the top panel of Figure 5 
after converting redshift to look-back time in a cosmology compatible with stellar ages. The lower 
panel in Figure 5 shows the age-metallicity relationship for disk stars determined in the landmark 
study by Edvardsson et al. (1993). This sample includes stars with the kinematics of both 
thin and thick disk, defined in terms of the mean velocity perpendicular to the Galactic plane: 
{\W\) = 19 km s~^ and ^ 42 km for thin and thick disk stars respectively (Freeman 1991). 
In constructing their sample, Edvardsson et al. aimed to include approximately equal numbers 
of stars in each metallicity bin above Z = O.IZq; consequently, metal-poor stars are relatively 
over-represented in Figure 5. 

Stellar ages are notoriously uncertain, as is the mapping of redshift to look-back time. 
However, even allowing for an arbitrary sliding of the points in Figure 5 along the x-axes, the 
metallicity measurements in DLAs evidently do not match the chemical evolution of the Milky 
Way disk. The typical value ^dla = —1.13 is lower than that of even the most metal-deficient 
stars in the Edvardsson et al. survey, and at all ages the spread of chemical abundances in the 
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disk is smaller than that of the DLA sample. 

This point is reinforced by Figure 6, where we compare the metallicity distribution of DLAs 
with those of stars in the thin disk, thick disk and halo populations; in Figure 7 we show the 
comparison with the metallicity histogram for globular clusters. Values for disk stars are from the 
work by Wyse & Gilmore (1995). These authors combined spectroscopic determinations of [Fe/H] 
for a sample of F and G stars located 1 — 5 kpc from the plane with data for samples near the Sun, 
paying particular attention to including only stars with potential main-sequence lifetimes greater 
than 12 Gyr. That is, their combined sample should be complete, in the sense of not missing 
disk stars which have by now evolved away from the main sequence, and the resulting metallicity 
distributions presumably provide an integrated record of the chemical evolution of the disk. The 
thin disk distribution in Figure 6 includes the low metallicity tail discussed by Wyse &: Gilmore 
(1995); similarly, the thick disk histogram is consistent with the metal- weak tail shown in Figure 
22 of Beers & Sommer-Larsen (1995). The halo sample is from the survey of high proper-motion 
stars in the solar neighbourhood by Laird et al. (1988), while the histogram in Figure 7 is based 
on the distribution of [Fe/H] in 40 globular clusters plotted in Figure 16 of Carney et al. (1996). 

The comparison between the metallicity distribution of DLAs and those of stellar populations 
in the Galaxy is complicated by the fact that about half of the values which make up the bins 
with [Zn/H]DLA < —0.8 in Figures 6 and 7 correspond to upper limits of [Zn/H] in our survey. 
Were we to exclude the upper limits from the sample, the resulting distribution would be skewed 
to higher metallicities. This is also the case if they are included in the sample as detections, as 
we have done; therefore the true distribution of Zdla is both broader and shifted towards lower 
metallicities (by undetermined amounts) than the histogram reproduced in Figures 6 and 7. 

Bearing this in mind, the middle and bottom panels of Figure 6 show that the metallicity 
distribution of DLA galaxies is different from those of long-lived stars in the Galactic disk. 
Although there is some overlap with the thick disk histogram, the bulk of stars in the disk of the 
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Milky Way apparently formed from gas which was significantly more metal-rich than that giving 
rise to damped Lyman a systems. The narrow distributions for disk stars in Figure 6 reflect the 
finding by Edvardsson et al. (1993) that the average metallicity has increased very little over the 
lifetime of the disk; the scatter at any age in the bottom panel of Figure 5 is nearly as large as the 
difference in mean metallicity over the entire time span considered. This is also the case for the 
old open clusters of the Milky Way disk (Friel 1995). 

The width of the ^dla distribution is comparable to those of halo stars and globular clusters, 
but it peaks at a higher metallicity. This is probably a real effect, rather than being due to the 
inclusion of upper limits in our sample (as discussed above), since the column density weighted 
mean metallicity is [(Zn/H^LA)] = —1.13. We consider it unlikely that the offset between the 
observed and true peaks of the Zdla distribution is as large as required to bring the histograms 
in the top panel of Figure 6 and in Figure 7 into agreement. Rather we favour the interpretation 
that, as a whole, the population of DLA galaxies is genuinely more metal enriched than the stellar 
components of the Galactic halo. 

The comparisons discussed above lead to two possible conclusions concerning the nature of 
damped Lyman a galaxies. The most straightforward interpretation is that a wide range of galaxy 
morphological types, at different stages of chemical evolution, make up the the DLA population. 
Available imaging data at 2; ^ 1 are certainly consistent with this view. A more intriguing 
possibility is that DLA systems at high redshift arise primarily in the spheroidal component of the 
present-day galaxy population, by analogy with the interpretation of the U drop-out galaxies put 
forward by Steidel et al. (1995, 1996). In the Milky Way, the halo and inner bulge may well be 
related, with the halo having lost 90% of its mass to the bulge (e.g. Wyse, Gilmore, h Franx 
1997); in this picture the halo-bulge system is an evolutionary sequence parallel to that of the 
thick disk-thin disk. One could speculate, then, that the distribution of .^dla, with its peak at a 
higher metallicity than halo stars and globular clusters, reflects different stages in the transition 
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from metal-poor halo to a predominantly metal-rich bulge (Ibata &; Gilmore 1995). 

4.2.1. Divergent Clues from the Absorption Line Profiles? 

The message conveyed by Figure 6 contrasts with the interpretation by Wolfe and collaborators 
of the complex absorption line profiles, often extending over more than 100 km s~^, revealed by 
high resolution spectroscopy of DLAs (Wolfe 1995; Prochaska & Wolfe 1997b). These authors have 
argued that in many cases the different components which make up the absorption lines are not 
distributed at random in velocity; rather, there appears to be a more regular trend of decreasing 
optical depth with increasing velocity difference from the wavelength where the absorption is 
strongest. This 'edge-leading asymmetry' is the pattern which would be produced by a rotating 
thick disk, intersected at some distance from the centre, if the average density of gas falls off 
with distance from the centre and from midplane. Prochaska 8z Wolfe show that the frequency 
with which such absorption profiles are encountered is consistent with expectations for randomly 
oriented disks; this leads them to conclude that most, if not all, DLAs arise in large {R > 10 kpc) 
disks with high rotation velocities {vrot ~ 200 km s""-*^). Such structures, if common at z ^ 2, are 
very difficult to explain in currently favoured models of galaxy formation (e.g. Baugh et al. 1997). 

The Milky Way is the only galaxy for which we have a record of both chemical abundances 
and kinematics over its past history. Based on this body of data, the metallicitics wc measure 
in the damped Lyman a systems appear incompatible with the rotating disk interpretation 
put forward by Prochaska & Wolfe. This can be appreciated by considering compilations of 
metallicitics and velocities now available for large samples of stars, such as that published recently 
by Carney et al. (1996). From their Figures 1 and 3 it can be seen that, of the stars with metal 
abundances similar to those of DLAs, approximately half have retrograde motions; at a metallicity 
Z = —1.1 the mean velocity relative to the disk rotation is (V) ~ —150 km s~^ . This point is best 
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illustrated by Figure 5 of Carney et al. which shows the metaUicity histograms in various intervals 
of V; our distribution of Zdla corresponds to values of V in the -100 to -200 km s'^. 

Evidently, when our Galaxy had an average metallicity of ^ 1/10 of solar, it did not exhibit the 
kinematics of a disk rotating at ~ 200 km s~^. 

Reconciling these contrasting clues to the nature of damped Lyman a galaxies is an important 
task for the future. Here we put forward three possible ways out of the current impasse: 

1. Our Galaxy is atypical, and the physical processes which gave rise to its stellar populations 
were not shared by most other galaxies at high redshifts. Although this possibility cannot be 
discounted, it is not a very constructive hypothesis to take refuge in, as it will be difficult to test 
it observationally — at least in the near future. 

2. The absorption profiles are being overinterpreted. A possible concern here is that material 
whose motion is due not to rotation but to energetic events, such as supernova shocks, may 
contribute to the ultraviolet absorption lines, since these transitions are sensitive to even relatively 
small column densities of gas. The 'edge-leading asymmetry' interpretation was first proposed 
by Lanzetta & Bowen (1992) in their analysis of 13 Mg II absorption components spread over 
250 km s~^ in the Zabs = 0.39498 DLA in Q1229— 021. However, it is far from clear that this is 
really a massive disk; from their analysis of iJST images of the field, Le Brun et al. (1997) propose 
that the absorber is instead a faint {Lb < O.IL*) low surface brightness galaxy. Furthermore, 
strong Mg II absorption spanning ^ 300 km s~^ can also be produced by galaxies which are nearly 
face-on, such as M61 (Bowen, Blades, & Pettini 1996). All these factors cast some doubts on 

a detailed correspondence between the profiles of ultraviolet absorption lines and the large-scale 
kinematics of the intervening galaxies. 

3. A third option, and one which we have already proposed, is that DLA galaxies comprise a 
mix of different morphological types. Thus, it is conceivable that some do exhibit the kinematics 
of rapidly rotating disks, while others may be spheroids or irregular star-forming galaxies with less 
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ordered velocity fields. This is a hypothesis which can be tested. As more cases become available 
where both kinematics and chemical abundances are measured in the same DLA, it will be of 
great interest to examine whether there is any correlation between these two parameters, as found 
in the stellar populations of the Milky Way. 

5. SUMMARY AND SUGGESTIONS FOR FUTURE WORK 

We have assembled the largest sample of damped Lyman a systems for which metallicities 
have been measured free from the complications introduced by dust depletions. The expanded 
data set reinforces the two main conclusions reached in our earlier study (Pettini et al. 1994): 
(1) DLAs are generally metal-poor, at all redshifts sampled; and (2) there is a large spread in 
abundances at all epochs. We interpret these findings as evidence for a protracted epoch of galaxy 
formation, and propose that galaxies of different morphological types and at different stages of 
chemical evolution make up the DLA population. 

The metallicity distribution of DLAs is broader and peaks at lower metallicities than those of 
either the thin or thick disk of our Galaxy. Thus, the chemical abundance data presented here do 
not support the interpretation of the absorption line profiles in terms of thick disks with rotation 
velocities Vrot ~ 200 km s~^ most recently discussed by Prochaska & Wolfe (1997b). This apparent 
discrepancy may be resolved by further work on both the kinematics and the abundances. With 
the near-infrared spectrographs now being built for 8-10 m telescopes it will be possible to detect 
the familiar optical emission lines from star-forming regions in the absorbing galaxies. The widths 
of these features are likely to be more representative of the global kinematics than the ultraviolet 
absorption lines which can be so easily affected by local phenomena such as interstellar shocks. 
On the abundance front, the ratios of chemical elements manufactured in different nucleosynthetic 
processes have been used to good effect in unravelling the history of star formation in our Galaxy; 
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the same techniques are now beginning to be applied to high redshift DLAs (Pettini et al. 1995b; 
Lu et al. 1996). 

The column density weighted mean metallicity of DLAs at 2; ^ 2 is in agreement with 
expectations based on the metal ejection rate deduced by Madau (1996) from the integrated 
ultraviolet luminosity of star forming galaxies at these redshifts. Our data, when combined with 
the [Fe/H] measurements by Lu et al. (1996), appear to reflect the rapid increase in the comoving 
star-formation rate between 2; 4 and ^ 2 indicated by the relative numbers of B and U drop-outs 
in the Hubble Deep Field. While these comparisons are of necessity still very approximate, the 
implication seems to be that observations of DLAs provide a reasonably accurate census of metal 
enrichment at these epochs. It is encouraging that three independent methods which have been 
applied to the quest for the epoch of galaxy formation — the global star formation rate deduced 
from the ultraviolet luminosity of high-redshift galaxies, the rate of consumption of neutral gas 
implied by the redshift evolution of J^dla, and the metallicity of DLAs — apparently give a broadly 
consistent picture of the universe at 2; ^ 2 . 

This is not the case at z ^ 1.5, where ^dla apparently does not rise as expected from simple 
models of cosmic chemical evolution. There are a number of plausible explanations for this, 
including the effects of dust, as discussed extensively by Fall and collaborators, and an increasing 
contribution of low surface brightness galaxies to the cross-section for DLA absorption. The major 
obstacle to progress in this area is still the paucity of DLAs with measured element abundances 
at intermediate redshifts. And yet it is essential to follow the evolution of the DLA population 
to the present time in order to be confident of our interpretation of the high redshift data. New 
DLAs at z ^ 1 are still being identified and the sample is slowly growing. With STIS on the 
HST measurements of [Zn/H] can be extended to redshifts lower than the limit Zabs — 0-65 of 
ground-based observations. In the next few years the 2dF and Sloan sky surveys (Taylor 1995; 
Gunn &: Weinberg 1995) are expected to increase the number of known DLAs by one order of 
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magnitude. With 8-10 m telescopes it will then be possible to repeat surveys such as this one 
towards substantially fainter, and potentially more reddened, QSOs. Such programmes should 
lead to a better assessment of the significance of dust bias in current DLA samples. Finally, with 
large telescopes we will soon be able to measure element abundances from the optical emission 
lines of galaxies at redshifts z ^ 0.1 — 0.5 . Such data will complement in a very important way 
the information provided by galaxies selected from their absorption cross-section. 
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TABLE 1 

DAMPED Lya SYSTEMS OBSERVED 



QSO 


V 




Ref. 




Telescope 


Resolution 


Integration 


S/N 






(mag) 










(A) 


Time (s) 




(mA) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0000-263 


18 


4.111 


1 


3.3901 


AAT 


1.1 


58 200 


30 


25 


0056+014 


18.9 


3.154 


2 


2.7771 


WHT 


1.0 


21000 


28 


28 


0836+113 


19.4 


2.696 


3 


2.4651 


WHT 


0.77 


34000 


14 


47 


0841+129 


17 


2.5: 


4 


2.3745 


WHT 


0.76 


20000 


30 


23 


0841+129 


17 


2.5: 


4 


2.4764 


WHT 


0.76 


20000 


30 


22 


0913+072 


17.1 


2.785 


5 


2.6183 


WHT 


0.75 


19800 


46 


14 


1151+068 


18.8 


2.762 


6 


1.7736 


WHT 


0.81 


14400 


44 


20 


1209+093 


18.5 


3.297 


7 


2.5843 


WHT 


0.82 


18000 


29 


24 


1337+113 


18.2 


2.919 


3 


2.7957 


WHT 


0.75 


18000 


9 


66 


2239-386 


18 


3.511 


8 


3.2810 


AAT 


1.1 


65 875 


18 


43 



3(T detection limit for the rest frame equivalent width of an unresolved absorption Hne. 

REFERENCES — (1) SavagHo, D'Odorico, & Moller 1994; (2) Wolfe et al. 1995; (3) Smith, Cohen, & Bradley 1986; (4) Hazard 1994; 
(5) Sargent, Boksenberg & Steidel 1988; (6) Turnshek et al. 1989; (7) Hazard, McMahon, & Morton 1987; (8) Lu & Wolfe 1994. 



TABLE 2 

REDSHIFTS AND EQUIVALENT WIDTHS OF Zn II AND Cr II LINES 



QSO 



(1) 



Line 1 
Zn II A2025.483 



(2) 



Wo 
(mA) 
(3) 



Line 2 
Cr II A2055.596 



(4) 



Wo 
(mA) 
(5) 



Line 3 
Cr II A2061.575+ 
Zn II A2062.005 

Wo 
(mA) 
(6) 



Line 4 
Cr II A2065.501 



(7) 



Wo 
(mA) 
(8) 



0000-263 
0056+014 
0841+129 
1151+068 
1209+093 



< 25 
1 1 

< 26 
1.7737 49 + 9 
2.5842 190+ 10 



3.3901 
2.7766: 
2.3746 
1.7736 
2.5846 



31 + 8 
116+17 
37+7 
78 + 8 
185+15 



22 + 7 
117+16 
52 + 7 
75 + 8 
195+15 



2.7769: 

1.7736 
2.5842 



< 25 
49+12 

< 23 
34 + 6 
120+15 



Notes. — Absorption redshifts are vacuum heliocentric. Equivalent widths are rest-frame values with +1(t errors. 
No redshifts are listed for Line 3 as two absorption lines contribute to the feature. 
^ Blended with atmospheric A band 



TABLE 3 

Zn AND Cr ABUNDANCES IN THE DAMPED Lya SYSTEMS 



QSO 

(1) 



(2) 



7V(Zn+) N{Zn+)/N(H°) Zn Underabun. iV(Cr+) N{Ct+)/N(E°) Cr Underabun. iV(Cr+)/7V(Zn+) 

(1012 cm-2) (10-9 ) Factor (lO^^ cm-^) (IQ-' ) Factor 

(3) (4) (5) (6) (7) (8) (9) 



0000-263 

0056+014 

0201+365" 

0302-223' 

0454+039'^ 

0836+113 

0841+129 

0841+129 

0913+072 

0935+417'' 

1104-180' 

1151+068 

1209+093 

1328+307^ 

1337+113 

1946+769" 

2239-386 



3.3901 

2.7771 

2.462 

1.0093 

0.8596 

2.4651 

2.3745 

2.4764 

2.6183 

1.3726 

1.6616 

1.7736 

2.5843 

0.6922 

2.7957 

2.8443 

3.2810 



<1.4 

3.5+2 

5.9+0.6 

3.0+0.6 

3.8+0.7 

< 2.6 
1.8+0.5 

< 1-2 

< 0.8 
1.8+0.4 
4.2+1.1 
2.5+0.5 
> 10 
5.3 

< 3.7 
<1.4 
<2.4 



<0.56 

2.7+1.6 

24.5+3 

14+4 

6.7+1.3 

< 7 

2.0+0.7 

< 2.1 

< 3.3 
7+2 
7+2 

1.25+0.4 
> 5 
2.8 

< 4.6 
<6.8 
<4.2 



>80 



1.85+0.2 



o 9+1.3 



> 6.5 

> 22 

> 14 



36ir 
< 9 
16 

> 10 
>6.7 
>11 



8+1.5 

28+4 

14.5+1 

13+2.5 

27+3 

< 12 
9.5+2 

< 5.6 

< 3.4 
15+4 
15+5 
22+2 
> 35 
20 

< 17 
<9.6 
<22 



3.1+0.9 

22+5 
60+7 
60+15 
47+5 

< 32 
10.5+3 

< 9.3 

< 15 
60+20 
25+9 
11+3 
> 17.5 
10.5 

< 21 
<48 
<38 



155^™ 



22^^ 



8.0+0.9 

ioi; = 

> 15 



45+19 



> 51 

> 32 



=+4 

'-2 



19^11 



4^+60 
^''-34 



< 27 
46 

> 23 
>10 
>13 



>5.5 

8+5 

2.5+0.3 

4.3+1.2 

7+1.5 

5.5+2.5 



8.5+3 

3.5+1.5 

9+2 

3.8 



CO 
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NOTES — 

" Prochaska k \No\ic 1996 

* Pettini & Bowcn 1997 

" Steidel et al. 1995 

^ Meyer, Lanzetta, & Wolfe 1995 

" Smette et al. 1995 

f Meyer & York 1992 

9 Lu et al. 1995 
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TABLE 4 

H I COLUMN DENSITIES IN THE DAMPED Lyo SYSTEMS 



QSO 2abs iV(H°) FWHM(LyQ) 

(1020 cm-2) (A) 
(1) (2) (3) (4) 



0000-263 

0056+014 

0201+365'' 

0302-223"^ 

0454+039"* 

0836+113 

0841+129 

0841+129 

0913+072 

0935+417" 

1104-180" 

1151+068 

1209+093 

1328+307'' 

1337+113 

1946+769' 

2239-386 



3.3901 

2.7771 

2.462 

1.0093 

0.8596 

2.4651 

2.3745 

2.4764 

2.6183 

1.3726 

1.6616 

1.7736 

2.5843 

0.6922 

2.7957 

2.8443 

3.2810 



25+5" 

13±2 

2.4+0.2 

2.15+0.35 

5.7+0.3 

3.8+0.4 

9+2 

6+1.5 

2.3+0.4 

2.5+0.5^ 

6 

20+5 

20+5 

19 

8+2 

2+0.5 

5.8-j 



18.0 

7.6 

10.6 
15.2 
12.5 
8.1 



22.0 
22.2 

14.0 



NOTES — 

" Savaglio, D'Odorico, & Moller 1994 

Prochaska & Wolfe 1996 
" Pettini & Bowen 1997 

Steidel et al. 1995 
" Meyer, Lanzetta, & Wolfe 1995 

Lanzetta & Meyer 1996, private communication 
« Smette et al. 1995 

Meyer & York 1992 
' Lu et al. 1995 
^ Lu & Wolfe 1994 
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Table 5. COLUMN DENSITY WEIGHTED METALLICITIES 





Redshift Range 


Lookback Time (Gyr)" 


DLAs 


Detections 


Upper Limits 


[(Z"/Hdla>] 


Pull Sample 


0.6922 - 


- 3.3901 


7.8 - 


- 14.7 


34 


19 


15 


-1.13 ±0.38 


Subsample 1 


0.50 - 


- 1.49 


6.3 - 


- 11.4 


4 


4 





-0.98 ± 0.33 


Subsample 2 


1.50 - 


- 1.99 


11.4 - 


- 12.7 


8 


6 


2 


-0.96 ± 0.44 


Subsample 3 


2.00 - 


' 2.49 


12.7 - 


- 13.6 


12 


6 


6 


-1.23 ±0.38 


Subsample 4 


2.50 - 


- 2.99 


13.6 - 


- 14.3 


7 


3 


4 


-1.11 ±0.27 


Subsample 5 


3.00 - 


- 3.49 


14.3 - 


- 14.8 


3 





3 


< -1.39 



''Ho = 50 km s-i Mpc-i; 50 = 0.01 
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Fig. la. — Portions of the QSO spectra observed in our survey of Zn II and Cr II lines in damped 
Lyman a systems. The x-axis is wavelength in A; the y-axis is residual intensity. The vertical tick 
marks indicate the expected positions of the absorption lines, whether they have been detected or 
not. Line 1 : Zn IT A2n25.48H: line 2: Cr IT A2n55.59fi: line .3: Cr IT A2nfil .575 + Zn IT A2nfi2.nn5 



-44- 



Q0913 + 072 z. =2.6183 



T 



1.0 
0.8 
0.6 



_L 



_L 



I 

2 



_L 



I I 
3 4 



7350 7400 7450 

Q1151+068 z. =1.7736 



T 



T 



1.0 
0.8 
0.6 




I 

2 



I 

3 



I 

4 



5650 5700 
Q1209 + 093 z^bs=2.5843 



—\ ' 1 '"^1 1 ' 1 ^ 



1.0 
0.8 
0.6 



I 

2 



I I 
3 4 

\ 



7250 



7300 7350 
Q1337+113 z^b,= 2.7957 



7400 




7700 



7750 7800 
Q2239-386 z,b,=3.2810 



7850 




0.6 - 



8650 



8700 



J 

8750 




8800 



8850 



-45- 



Q0056+014 z^^=Z.lll\ 




4300 4400 4500 4600 4700 4800 
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Q0841+129 z. =2.3745, 2.4764 
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3200 3300 3400 3500 3600 3700 



Fig. 2a. — Normalised portions of the blue spectra of QSOs in our survey centred on the damped 
Lyman a absorption line. The x-axis is wavelength in A; the y-axis is residual intensity. In each 
Danel the short-dash line shows the theoretical damoine orofile corresDondine to the value of neutral 
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Q1209+093 z,b.=2.5843 
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Fig. 2b . — (continued) 
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Fig. 3. — The abundance of Zn for the 34 damped Lyman a systems m the present survey pfotted 
against redshift. Abundances are measured on a log scale relative to the solar value shown by the 
broken line at [Zn/H] = 0.0. Upper limits, corresponding to non-detection of the Zn II lines, are 
indicated by downward-pointing arrows. For two damped systems, indicated by dots with upward- 
pointing arrows, we derive lower limits to the abundances, because the absorption lines may be 
saturated. 
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Fig. 4. — Cosmic chemical evolution as deduced from the metallicity of damped Lyman a systems 
(filled squares) and from the metal ejection rate, p'z, of high redshift galaxies in the Hubble Deep 
Field (dots). The filled squares are the values of [(Zn/Hj-,L/j^)] from the present survey; the broken 
line shows the predictions of a simple integration of the redshift evolution of p'z from Madau (1996) 
for a Hq = 50 km Mpc~^ and qo = 0.5 cosmology. For qq = 0.01 the broken line would shift to 
lower values of log {Z)/Zq by a factor of ~ 2 . 
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Fig. 5. — Top Panel: Available measurements of metallicity in damped Lyman a systems plotted 
as a function of lookback time for Hq = 50 km s~^ Mpc and qo = 0.01 . Bottom Panel: Metallicities 
of 182 F and G dwarf stars in the Galactic disk with measured iron abundances and ages from the 
large-scale study by Edvardsson et al. (1993). 
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Fig. 6. — Metallicity distributions, normalised to unity, of damped Lyman a systems and of stars 
belonging to the disk and halo populations in the Milky Way. See text for references to the sources 
of stellar data. Upper limits to [Zu/HJdla have been included as measurements in the histogram. 
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Fig. 7. — Metallicity distributions, normalised to unity, of damped Lyman a systems and of 40 
globular clusters in the sample of Carney et al. (1996). Upper limits to [Zn/H]DLA have been 
included as measurements in the histogram. 



